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Abstract 

We investigate singlet scalar dark matter (DM) candidate in a left-right symmetric gauge model with two 
Higgs bidoublets (2HBDM) in which the stabilization of the DM particle is induced by the discrete symme- 
tries P and CP. According to the observed DM abundance, we predict the DM direct and indirect detection 
cross sections for the DM mass range from 10 GeV to 500 GeV. We show that the DM indirect detection 
cross section is not sensitive to the light Higgs mixing and Yukawa couplings except the resonance regions. 
The predicted spin-independent DM-nucleon elastic scattering cross section is found to be significantly de- 
pendent on the above two factors. Our results show that the future DM direct search experiments can cover 
the most parts of the allowed parameter space. The PAMELA antiproton data can only exclude two very 
narrow regions in the 2HBDM. It is very difficult to detect the DM direct or indirect signals in the resonance 
regions due to the Breit-Wigner resonance effect. 
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I. INTRODUCTION 



The existence of dark matter (DM) is by now well established from astrophysical observations 
ti|]. Together with the recent WMAP results, the cosmological observations have shown that the 
present Universe consists of about 73% dark energy, 23% dark matter, and 4% baryonic matter yfl. 
In the standard model (SM) of particle physics, there is no cold DM candidate. Therefore, one has 
to extend the SM to account for the existence of DM. The DM candidate is often accompanied 
by some discrete symmetries to keep it stable, such as the R parity in super symmetric (SUSY) 
models and KK parity in universal extra dimension models. Although the discrete symmetries are 
necessary for the DM stability, they may be introduced from different motivations yj]. 

In the left-right (LR) symmetric gauge model |3|-l5|] with spontaneous CP violation (SCPV), the 
P and CP symmetries are exact before the spontaneous symmetry breaking (SSB). In this case, it 
is possible that the discrete symmetries P and CP strongly constrain the scalar sector of the model 
and naturally give stable DM candidates. This possibility has not been emphasized in the literature, 
due to the fact that most of the popular models such as SM and SUSY violate P maximally. In 
Ref. [6], we have shown that the P and CP symmetries can give a stable DM candidate in an 
extension of a left-right symmetric gauge model with a singlet scalar field S = (S^ + iS D )/ V2. In 
this model, the CP odd particle S D is stable even after the SSB, provided that it does not develop 
vacuum expectation value (VEV). 

Without large fine-tuning, it is difficult to have a successful SCPV in the minimal left-right 
symmetric gauge model with only one Higgs bidoublet (1HBDM) fl |7|]. This is because in the 
decoupling limit the predicted CP violating quantity sin 2/3 ~ 0.1 with ft being a CP phase angle in 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix is far below the experimentally measured value 
of sin 2/3 = 0.671 + 0.024 from the two B-factories [8]. In addition, the 1HBDM is also subject 
to strong phenomenological constraints from low energy flavor changing neutral current (FCNC) 
processes, especially the neutral kaon mixing which pushes the masses of the right-handed gauge 
bosons and some neutral Higgs bosons much above the TeV scale Motivated by the require- 
ment of both spontaneous P and CP violations, we have considered the left-right symmetric gauge 
model with two Higgs bidoublets (2HBDM) [10]. In the 2HBDM, the additional Higgs bidoublet 
modifies the Higgs potential so that the fine-tuning problem in the SCPV can be avoided, and the 
bounds from the FCNC processes can be relaxed. The extra Higgs bidoublet may also change the 
interferences among different contributions in the neutral meson mixings, and lower the bounds 
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for the right-handed gauge boson masses not to be much higher than the TeV scale [10]. Such a 
right-handed gauge boson can be searched at the LHC using the angular distributions of top quarks 
and the leptons from top quark decays lllin . 

In Ref. B6D, we have shown that the discrete symmetries P and CP can be used to stabilize 
the DM candidate S D in the 1HBDM and 2HBDM with the SCPV. Using the observed DM abun- 
dance, we can constrain the parameter space and predict the spin-independent (SI) DM-nucleon 
elastic scattering cross section. For simplicity, we have only considered the case with no mixing 
among light neutral Higgs bosons in the 2HBDM and the dark matter is heavy. In this paper, we 
shall demonstrate in detail the mixing effect on the DM direct detection. Notice that several new 
DM annihilation channels can be derived, namely two DM particles may annihilate into a gauge 
boson and a Higgs boson. On the other hand, we are going to extend the DM mass range from 
200 GeV < m D < 500 GeV to 10 GeV < mo < 500 GeV. As a consequence, one will meet several 
resonances in the 2HBDM. Therefore we shall consider the Breit-Wigner resonance effect for the 
determination of the DM relic density uM- In addition, we will also consider the DM indirect 
search in the 1HBDM and 2HBDM. The paper is organized as follows: In Section. HO we outline 
the main features of the 1HBDM and 2HBDM with a singlet scalar. In Sec. [Ill] and Sec. [IV] we 
discuss the parameter space, the DM direct search and the DM indirect search in the 1HBDM and 
2HBDM, respectively. Some conclusions are given in Sec. |V] 



II. THE LEFT-RIGHT SYMMETRIC GAUGE MODEL WITH A SINGLET SCALAR 



We begin with a brief review of the 2HBDM described in Ref. Q]. The model is a simple 
extension to the 1HBDM, which is based on the gauge group SU(2) L <8> SU(2) R <g> U(\) B -l- The 
left- and right-handed fermions belong to SU(2) L and S U(2) R doublets, respectively. The Higgs 
sector contains two Higgs bidoublets <p (2,2* ,0), x (2,2* ,0) and a left (right) -handed Higgs triplet 



Ar,(R) (3(1),1(3),2) with the following flavor contents 
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The introduction of Higgs bidoublets <p and^- can account for the electroweak symmetry breaking 
and overcome the fine-tuning problem in generating the SCPV in the 1HBDM. Meanwhile it also 
relaxes the severe low energy phenomenological constraints lioll . Motivated by the spontaneous P 
and CP violations, we require P and CP invariance of the Lagrangian, which strongly restricts the 



structure of the Higgs potential. The most general potential containing only the (p and A LiR fields 
is given by 

% A = -^Tr(0V) - ^[Tr(0V) + Tr(^)] - ^[Tr(A L A[) + Tr(A R A^)] 
+A 1 m^<t>)] 2 + A 2 {[Tr(^cf>)] 2 + [Tr(^)] 2 } + A 3 [Tr(fcf>)Tr(^)] 
+/t 4 {Tr(0V)[Tr(0V) + Tr^)]} 

+ Pl {[Tr(A L A[)] 2 + [Tr(A R A^)] 2 } + p 2 [Tr(A L A L )Tr(A[A[) + Tr(A R A R )Tr(A^)] 
+p 3 [Tr(A L A[)Tr(A R A;)] +p 4 [Tr(A L A L )Tr(A;A^) + Tr(A[A[)Tr(A R A R )] 
+a 1 Tr(0V)[Tr(A L A[) + Tr(A R A^)] + a 2 Tr[(0V) + (# + )]Tr[(A L A[) + (A R A£)] 
+a 3 [Tr(00 t A L A[) + Tr(<^0A R A^)] 

+ y 8 1 [Tr(0A R t A[) + Tr^A^A^)] + /? 2 [Tr(0A R t A[) + Tr(^A L ^)] 
+A[Tr(0A R t A[) + Tr(0^A L 0Ag)], (2) 

where the coefficients Aj, p,, a, and in the potential are all real as all the terms are self- 
Hermitian. The Higgs potential "V^a involving ;f field can be obtained by the replacement^ <-» cf> 
in Eq. ©. The mixing term 'V^a can be obtained by replacing one of 4> by x in all the possible 
ways in Eq. (O. In order to simplify the discussion, we shall first consider the 1HBDM which 
already contains the main features of the complete model. Then we postpone the discussions on 
the^ contributions into Section ITVl 

After the SSB, the Higgs multiplets obtain nonzero VEVs 

« 2 > = % and (8%) = ^ , (3) 

where K\, /c 2 , vl and vr are in general complex, and k = Vkil 2 + |k 2 | 2 « 246 GeV represents the 
electroweak symmetry breaking scale. Due to the freedom of gauge symmetry transformation, 
one can take K\ and v R to be real. To avoid the fine-tuning problem of fermion masses, we require 
v L - and k 2 <K K\. The value of v R sets the scale of left-right symmetry breaking which is directly 
linked to the right-handed gauge boson masses, vr is subjected to strong constraints from the K, 



B meson mixings [|4j,l8|, 



9] 



14fl . The kaon mass 



as well as low energy electroweak interactions H13L 
difference and the indirect CP violation quantity 6k set a bound for vr around 10 TeV [|13L 

In our model, the P and CP symmetries have been required to be exactly conserved before the 
SSB, thus the discrete symmetries P and CP can be used to stabilize the DM candidate. In the 
framework of 2HBDM with a complex singlet scalar S = (So- + iS D )/ V2, we have considered 





p 


CP 




P 


CP 




P 


CP 




ft 


<t>* 


s +s* 


+ 


+ 


s -s* 


+ 




X 




X* 


ss* 


+ 


+ 


Tr(^<f>) 


+ 


+ 


A L(R) 




A* 

L(R) 




+ 


+ 


Tr((f>^(f) - ft<p) 






s 


S 


s* 


Tr(A[A L + A^A«) 


+ 


+ 


Tr(A[A L - A^Afi) 




+ 



TABLE I: The P and CP transformation properties of the Higgs particles and their gauge-invariant combi- 
nations. The "+" and "-" denote even and odd, respectively. 

this possibility in Ref. |(|]. The P and CP transformation properties of the Higgs particles and 
their gauge-invariant combinations have been shown in Table HI It is clear that the odd powers of 
(S - S*) are forbidden by the P and CP symmetries. Therefore So is a stable particle and can 
be the DM candidate when the VEV Va-f V2 of S is real. Although P and CP are both broken 
after the SSB, there is a CP type Z 2 discrete symmetry on S D remaining in the singlet sector. 
This discrete symmetry is induced from the original CP symmetry. We have checked that the 
P and CP transformation rules for S defined in Table U is actually the only possible way for the 
implementation of the DM candidate. 

For the annihilation cross section of approximately weak strength, we expect that the DM mass 
is in the range of a few GeV and a few hundred GeV. However, the mass m D of S d is related to the 
LR symmetry breaking scale vr ~ 10 TeV. To have a possible light DM mass, we may consider 
an approximate global U(l) symmetry on 5, i.e. S — » e l6 S . Then the P and CP invariant Higgs 
potential involving the singlet S is given by 



<V S = -fi 2 D SS* +A D (SS*) 2 + Y J A UD SS*O i -'^(S -S*) 2 , 



(4) 



where O x = Tr(0VX 2 = TrC^ + and 3 = Tr(A[ A L + A^A S ). Only the last term explicitly 
violates U(l) symmetry. After the SSB, S obtains a real VEV v a l V2. Then one can straightly 
derive 

% = ^[(Sl + 2v £r 5 cr + Si ) 2 - v 4 J + J] ^f( s l + 2V - S <r + + S 2 D )(.Oi - <0,» + ^S 2 D , (5) 



f=l 



where we have used the minimization condition fj. 2 D = A D vl. + from the singlet 5^ to 

eliminate the parameter The terms proportional to odd powers of So are absent in Eq. © 
which implies S D can only be produced by pairs. Notice that the mass term of S D should be 



absent with an exact global U(l) symmetry. As discussed in Ref. \^\, the explicit breaking of this 
U(l) symmetry can explain the naturalness of a light DM mass m D , but it does not destroy the 
stability of the DM candidate S D . 
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TABLE II: The mass spectrum for the Higgs and gauge bosons in the left-right symmetric gauge model with 
one Higgs bidoublet in the limit vl - and k% K\ . <f>. and <f>9' stand for real and imaginary components 
of (pi = (0? r + i(f>l')l V2, respectively. The gauge boson Z\{W\) corresponds to the Z(W) boson in the SM. 

The terms Ivo-So-Oj in Eq. © indicate that S a will mix with the Higgs bosons (p ( \ r , <p° 2 r , 6° L r and 
6° R r . The relevant mass matrix elements are given by 

M 2 = 2A D vl , M 2 = A hD KVa- , M 2 = 2A 2 , d kv ct , M 2 ^,. = A 3 ^ D v CT v L , M 2 ^,. = A^v^Vr . (6) 



For simplicity here we require v a > vr ~ 10 TeV » k which means the mixing angles between So- 
and the above four neutral Higgs bosons are small. The terms v 2 0, in Eq. © do not change the 
minimization condition forms for cp and A L(R ). This is because these terms only change the overall 
coefficients fi\, pi 2 and ^ in Eq. ©. Hence the mass matrixes of the Higgs multiplets (p and A LR 



160, which also indicates that the additional 



remain the same as that in the 1HBDM in Refs. J5, 
potential term <V S in Eq. © does not help in resolving the fine-tuning problem. Due to v L - 
and k 2 <sc K\, the mass eigenstates for the Higgs bidoublet and triplets approximately coincide with 
the corresponding flavor eigenstates. The mass spectrum for the Higgs and gauge bosons is listed 
in Table HI There is only one light SM-like Higgs h° from the real part of cp° r The masses of all 
the other scalars are set by vr which can be very heavy. From the Lagrangian in Eq. © one can 
easily obtain the interaction terms among the scalars. Some of the relevant cubic and quartic scalar 
interaction vertexes are listed in Table Hill 
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TABLE III: The cubic and quartic scalar vertexes among Higgs singlets and multiplets, where 
HH* stands for any states of (fc°fc ,fl°fl°,A°A°,flJfl£) and AA* stands for any states of 
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III. DARK MATTER SIGNAL IN THE 1HBDM 




FIG. 1: Feynman diagrams for the DM annihilation in the 1HBDM. 

As discussed in Sec. HB an approximate global U(l) symmetry on S can naturally lead to a 
light DM mass m D . Here we focus on 10 GeV < m D < 500 GeV. Considering the case v a > vr ~ 
10 TeV » k, one may find that most of scalar bosons in TablelUare very heavy except the SM-like 
one h°. In this case, the possible annihilation products are h°h°, W{W\IZiZ\ and fermion pairs // 
as shown in Fig. \T\ For ^-channel annihilation processes, the intermediate particles may be h°, H^, 

and H®. Because of v L - 0, one may neglect the case. In addition, the H° { contribution is 
also negligible as m H o » m h o. For the // annihilation process, the main contribution comes from 
the h° exchange diagram. This is because H% dominantly couples to the very heavy right-handed 
Majorana neutrinos (the corresponding annihilation process is kinematically forbidden). For the 
W\Wi/ZiZ\ processes, the diagram involving H® is suppressed by m H a » m h o. Notice that S a may 
be the intermediate particle for the h°h° case. It is clear that the dominant annihilation processes 
in Fig. [Dare the same as that in the minimal extension of SM with a real gauge singlet scalar when 



mo < m h o llnll . In the 1HBDM, the DM annihilation cross sections & = AE\E 2 crv (E\ and E 2 are 



the energies of two incoming DM particles) for different annihilation channels have the following 



forms: 
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where s is the squared center-of-mass energy [18]. The quantity F is defined as F(ij h o) 
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arctanh(^o)/^ ft o with g h o = ^js - 4m 2 ^Js - 4m 2 h0 /(s - 2m 2 ). The Higgs decay width T h o and 
G\ are given by 
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From Eqs. (1711111) seven unknown parameters enter the expression of total annihilation cross sec- 
tion, namely, m h o, m D , Ai D , a\Aj, D , m 2 , m\ and v R . For the mass of SM-like Higgs, we take 

2 

m h o = 120 GeV in the following parts. In fact, one may neglect the squared center-of-mass en- 
ergy s in the terms s - m 2 and s - in 2 , since the masses of s a and //? are around Vr. In a good 

2 

approximation, we find that only three independent parameters 



m D ,A hD and A R = ai/t 3 V( 2 Pi) 



(12) 



are relevant to our numerical analysis. Here we have used m 2 = 2p\V 2 R as it is shown in TableHO 



A. Constraints from the DM relic density 



In order to obtain the correct DM abundance, one should resolve the following Boltzmann 
equation 



(13) 
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where Y = n/s(x) denotes the DM number density. The entropy density s(x) and the Hubble 
parameter H evaluated at x = 1 are given by 

\4n>g* m 2 D 



s(x) = 



45 



45 M PL ' 



(14) 



where M PL 1 .22 x 10 19 GeV is the Planck energy, g t is the total number of effectively relativistic 
degrees of freedom. The numerical results of g t have been presented in Ref. 112011 . Here we take the 
QCD phase transition temperature to be 150 MeV. The thermal average of the annihilation cross 
section times the relative velocity (crv) is a key quantity in the determination of the DM cosmic 
relic abundance. We adopt the usual single-integral formula for (crv) 112 ill : 



1 mn r°° ,- xtJs 

(<tv) = — — J- o-( S )V?^i(— )d 
nl n 64n 4 x J 4m 2 m D 



n 2 64n 4 x 



(15) 



with 
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(16) 



2n 2 x '" \ s 

where K x {x) and K 2 (x) are the modified Bessel functions, x = m D /T and g, = 1 is the internal 
degrees of freedom for the scalar dark matter S d- In terms of the annihilation cross section & in 
Eqs. (TTIfTOl). one can numerically calculate the thermally averaged annihilation cross section (crv). 
Finally, we may obtain the DM relic density Q. D h 2 = 2.74 x 10 8 Y m D /GeV by use of the result 
Yo of the integration of Eq. (PT3l) . 

When the DM mass m D is larger than the mass of top quark, one will not meet the resonance 



III2I1 and threshold I122h effects in our model. Thus we use the approximate formulas to calculate the 



DM relic density for 200 GeV < mo < 500 GeV. In this case, (crv) can be expanded in powers of 
relative velocity and x~ x for nonrelativistic gases. To the first order (crv) - crox~ n , where n = 0(1) 
for s(p)-wave annihilation process fl 1 9TI . The approximate formula for (crv) is given by [12311 



(O-V) = CT X " = \ 



to - -C2co - co')x 1 + . . . 
2 



(17) 



where co = (d-ff + d- ZlZl + & WlWl + (T A o ft o)/4 and the prime denotes derivative with respect to s/4m 2 D . 
co and its derivative are all to be evaluated at s/4m 2 D = 1. Then Q. D h 2 is given by [19] 



(n + l)x n f +l 
a D h 2 = 1.07 x 10 9 — GeV" 1 

gj M PL cr 



(18) 



with 



x f = ln[0.038(n + l)(g i /g! /2 )M PL m D( r ] - (n + 1/2) ln{ln[0.038(n + l)(g,/g! /2 )M PL m D 6r ]}. (19) 




FIG. 2: Left panels: the predicted coupling A\o as a function of Ar and the DM mass nip from the observed 
DM abundance in the 1HBDM. Right panels: the predicted DM-nucleon scattering cross section cr^ 7 in the 
1HBDM with current and future experimental upper bounds. 

Notice that we take g, = 345/4 for 200 GeV < m D < 500 GeV. 

In terms of the observed DM abundance 0.1088 < Q. D h 2 < 0.1 158 B2J], we numerically solve the 
Boltzmann equation and derive the coupling Ai jD with different A R for 10 GeV < m D < 200 GeV. 
The numerical results are shown in Fig. [2] (upper- left panel). Due to the resonance contribution, 
a very small value of the coupling A\ D can be derived from the observed DM abundance for 
the resonance region (0.8 m h o < 2m D < m h o). Except for the resonance region, one may find 
Ai d ~ 0(10 2 - 10 _1 ). The parameter A R plays an important role to determine the DM relic 
density if m D > m h o. For illustration, we also plot the Ar = +0.1 cases which can significantly 
change the predicted A 1D as shown in Fig. [2] In fact, A\ tD may be very small (even to be zero) 
for the larger \A R \. In this case, the //^-exchange annihilation process is dominant. Here we have 
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assumed A\ t o is positive. If we simultaneously change the signs of A\, D and A R , the negative A\, D 
case may be approximately induced from the positive case. This feature can be well understood 
from Eqs. (1 1 OH 1 II) . It should be mentioned that the thermally averaged annihilation cross section 
(<xv) will significantly change as the evolution of the Universe when the DM particle is nearly one- 



half the mass of a resonance 111 211 . This is the Breit-Wigner resonance effect which has been used to 
explain the recent PAMELA [24J], ATIC Q25Q and Fermi Q26D anomalies. Notice that the decaying 
S d with a lifetime around O(\0 26 s) can also account for the electron and positron anomalies Il27ll . 
Here we have considered the Breit-Wigner resonance effect for the determination of the coupling 

For 200 GeV < m D < 500 GeV, we use the approximate formulas to scan the whole parameter 
space A\ D and A R . The allowed parameter space is shown in Fig. [2] (lower-left panel), which gives 
an allowed range -0.17 < A\ D < 0.17 and -0.32 < Ar < 0.32. The central region of this figure 
is excluded since these points can not provide large enough annihilation cross section to give the 
desired DM abundance. Notice that the approximate global symmetry U(l) requires rn^/vj, <sc A\^ 
which means the region near A\^ D = is disfavored. 

B. Dark matter direct search 

For the scalar dark matter, the DM elastic scattering cross section on a nucleon is spin- 
independent, which is given by 

SI 4/ m D m n \ 2 (Zf p + (A-Z)f n ) 2 



v n I I ~ \ 2 , (20) 

n\m D + m n ) A 2 

where m n is the nucleon mass. Z and A -Z are the numbers of protons and neutrons in the nucleus. 
f p n is the coupling between DM and protons or neutrons, given by 

fp,n = J] f ?T a <lit + ^ ^"^7' 

q=u,d,s 1 q=c,b,t 1 

where ff u = 0.020 + 0.004, ff d = 0.026 + 0.005, = 0.118 + 0.062, fj£ = 0.014 + 0.003, 
fi] = 0.036 + 0.008 and = 0.118 + 0.062 J28]. The coupling ffy n) between DM and gluons 
from heavy quark loops is obtained from f^y = 1 - 2 9 =»,d,.v/^' l) 5 which leads to ~ 0.84 and 
f% * 0.83. In our model, the DM-quark coupling a q in Eq. (|2T|) is given by 



X\n m 



a q = -^-f • (22) 
2m D m z h0 
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Because of f n w f p , we can derive 



r 



(23) 



7T \m£» + m n/ 

It is worthwhile to stress that erf 7 is independent of Ar. 

Using the predicted A\p from the observed DM abundance, we straightly calculate the spin- 
independent DM-nucleon elastic scattering cross section <xf 7 . The numerical results are shown in 
Fig. [2] (right panels). For lOGeV < m D < 200 GeV, we find that two DM mass ranges can be 
excluded by the current DM direct detection experiments CDMS II [|29D and XENON 10 [301]. Due 
to the existence of Ar, we can obtain different values of erf/ for a given DM mass m D when the 
annihilation channel S dS d —* h°h° is open. In this case, one can obtain erf/ < 7 x 10~ 45 cm 2 for 
200 GeV < m D < 500 GeV as shown in Fig. [2] (lower-right panel), which is below the current 
experimental upper bounds. Nevertheless the future experiments XENON100 113 ill . CDMS 100 
kg 113 211 and XENON1T (13311 can cover most parts of the allowed parameter space. For the region 
near the resonance point, the predicted erf 1 is far below the current and future experimental upper 
bounds. 



C. Dark matter indirect search 

As shown in Sec. IIII Al (crv) is a key quantity in the determination of the DM cosmic relic 
abundance. On the other hand, (crv) also determines the DM annihilation rate in the galactic 
halo. It should be mentioned that the DM annihilation in the galactic halo occurs at v ~ 10~ 3 
(x « 3/v 2 = 3 x 10 6 ). Thus we calculate the thermally averaged annihilation cross section at 
x « 3 x 10 6 , namely (ctv)q. The numerical results have been shown in Fig. [3] for 10 GeV < 
m D < 200 GeV. Notice that we can derive the similar results for different values of A R . One 
may find 1 x 10" 26 cm 3 sec^ 1 < (ctv)q < 3 x 10~ 26 cm 3 sec~ J for most parts of the parameter 
space. The enhanced and suppressed (<xv)o on the two sides of the resonance point originate 
from the Breit-Wigner resonance effect lll2ll . When m D is slightly less than the W\ boson mass, 
the channel S £>S d — > W^WY 1S open at high temperature, which dominates the total thermally 
averaged annihilation cross section and determines the DM relic density. However this channel 
is forbidden in the galactic halo. Thus the threshold effect leads to a dip around W\ threshold 
[12211. When 200 GeV < m D < 500 GeV, one can obtain <crv) * 2.3 x 10^ 26 cm 3 sec 1 which 



is consistent with the usual s-wave annihilation cross section (crv) ~ 3 x 10 cm sec at the 
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freeze-out temperature Xf ~ 20. 
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FIG. 3: The predicted thermally averaged DM annihilation cross section (crv)o in the 1HBDM. 



In our model, the DM annihilation can generate primary antiprotons which can be detected 
by the DM indirect search experiments. Recently, the PAMELA collaboration reports that the 
observed antiproton data is consistent with the usual estimation value of the secondary antiproton 
||24f | . Therefore one can use the PAMELA antiproton measurements to constrain (crv)o- In Fig. 
[31 we have also shown the maximum allowed (crv)o for the MIN, MED and MAX antiproton 
propagation models given in Ref. Il34n . Then we can find that a very narrow region can be excluded 
by the PAMELA antiproton data in our model. In fact, the width of this excluded region is about 
0.4 GeV for the MED and MAX cases. When double DM mass 2m D is slightly less than the Higgs 
mass m h o, the predicted erf 1 and (<xv)o are very small which means that it is very difficult to detect 
the DM signals. 



IV. DARK MATTER SIGNAL IN THE 2HBDM 



We have discussed the Higgs singlet S D as the cold DM candidate in the 1HBDM. In this sec- 
tion, we generalize the previous discussions to the 2HBDM in which the other bidoublet y mixes 

n 

significantly with (p and A Lfi . In this case the SCPV can be easily realized HlOll . Comparing with 
the previous case, the main differences are that there could be more scalar particles entering the 
DM annihilation and scattering processes. Furthermore, the new contributions from these particles 
may modify the correlation between the DM annihilation and DM-nucleon elastic scattering cross 
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sections, which leads to significantly different predictions from the other singlet scalar DM models 
and the previous discussions. 



As shown in Eq. (OQ), the second Higgs bidoublet x contains two neutral Higgs contents x 



1,2" 



After the SSB, x\ 2 mav obtain the VEVs Wi >2 / V2. The squared sum of all the VEVs including 
k-1,2 should still lead to v E w = VVil 2 + N 2 + l^il 2 + \ w i? ~ 246 GeV. In general, the 2HBDM 
includes three light neutral Higgs bosons and a pair of charged light Higgs particles, whose masses 
are order of the electroweak energy scale. For simplicity, we consider k 2 ~ w 2 ~ 0. In this case, it 
is convenient for us to rotate Higgs bidoublets (p and^ into 



X = 



h 2 +ili3 r + 
X 2 



V2 



A) 



(24) 



o r 2 ) 

where are a pair of light charged Higgs bosons. Then one can diagonalize the mass matrix of 



three light neutral Higgs h\ 2 ^ and derive three light neutral Higgs mass eigenstates. The relation 
between h\ >2>3 and three mass eigenstates can be written as 



hi 
h 2 



c r c. 



s x c z 



-c x s y s z 



y -c x c y s z + s x s y -s x c y s z - c x s y c y c z 





( \ 




h 




H 




A 


) \ ) 



(25) 



where s x = sin 9 X , c x = cos 9 X and so on. Due to many unknown parameters in the Higgs potential 
of 2HBDM, we can not explicitly calculate three mixing angles 6 X , 8 y and 6 Z . For illustration, we 
consider three representative cases: (I) G x = 60°, 9 y = 60° and 9 Z = 150°; (II) 9 X = 30°, 9 y = 0° and 
9 Z = 0°; (III) 9 X = 0°, 9 y = 90° and 9 Z = 75°. The Case I means that there is the significant mixing 
among three light neutral Higgs. If all CP violation phases are absent, we can obtain 9 = 0° and 
9, = 0°. In the Case II, the light Higgs A is CP odd which does not mix with h and H. For the 
Case III, we only consider the scalar and pseudoscalar mixing, namely 9 X = 0°. 

In the 2HBDM, the possible DM annihilation products are //, W x W x IZ{Z X , WfH T \Z X (h, H, A), 
H + H~ and any two of the three neutral states (h,H,A) as shown in Fig. HI For a concrete nu- 
merical illustration, we choose all the masses m H , m A , m H ± = 180 GeV and m/, = 120 GeV. For 
cubic and quartic scalar vertexes, we assume they are the same as that in the 1HBDM. Namely, 
the vertexes of S D S D (h,H,A) and S D S D (h,H,A/H + )(h,H,A/H~) are set equal to -iAi D v E w and 
-iA\ D , respectively. Similarly, the cubic scalar vertexes among the light Higgs particles h, H,A 
and are set equal to -i3m 2 h /v E w, and the cubic scalar vertexes between S a and two light Higgs 
particles are assumed to be -iAi^Vo-. It is worthwhile to stress that the heavy Higgs particles from 
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So 



h,H,A/h+ S D 



h,H,A/h+ S D _ 



> -X 

/ h,H,A,H%,S a x 



So 



^h.H.A/h- S D - 



FIG. 4: Feynman diagrams for the DM annihilation in the 2HBDM. 

X' may be as the intermediate particles when two DM candidates annihilate into two light Higgs 
bosons. Nevertheless we still can use a coupling A R to describe the contributions of all possible 
heavy Higgs bosons. All annihilation cross sections & have been presented in Appendix lAl 
In the basis of Eq. ((24)) . the Yukawa interactions for quarks are given by 

- Ly = Q~l(yW + Y*& + *V + Y X X') Qr + h.c, (26) 

where Q L R = (u L j{, d L<R ) T . When both P and CP are required to be broken down spontaneously, 
the Yukawa coupling matrices Y*, Y^, Y x and Y x are complex symmetric. Then one may rotate the 
quark fields and derive the following Yukawa interactions relevant to light neutral Higgs particles: 
_ = h ll v^ L , v , + - , v + h_ ll p_ v , + h 1 -ih 1 - , + 
VI V / V2 V2 

where Y^' and Y$ are diagonal matrixes. According to the up and down quark masses, we can 
obtain Y^ q = V2m^/v EW and Y$ q = V2m 9 /v EW , respectively. In order to avoid the FCNC pro- 
cesses, we assume Y x ' and Y x ' are approximate diagonal matrixes due to approximate U(l) family 



symmetries [35] and require 

Y x m = R q Yf q and Y x qq = R q Y^ . (28) 

Since Y x ' and Y x ' don't contribute the quark masses, the parameter R q may be very large except 
the top quark case. 

In the 2HBDM, the parameter R q in Eq. (1281) controls the Yukawa couplings Y qq and Y qq . Fur- 
thermore, the parameter R q will affect the total annihilation cross section and change the predicted 
coupling A\ D . For illustration, we choose the following two scenarios 

R q = R = 1 and R q = R = 5 (q ± t and R, = 1 for the top quark) (29) 
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FIG. 5: The predicted coupling Ai t o and DM-nucleon scattering cross section cr n for three mixing cases 
in the 2HBDM with R = 1 and R = 5. 

to calculate the allowed coupling A\ D from the observed DM abundance. Considering three kinds 
of mixing cases and two R scenarios, we plot the allowed coupling A\ d for 10 GeV < m D < 
200 GeV in Fig. [5](left panels). It is clear that A\^ is dependent on the light Higgs mixing and the 
parameter R if m D < 120 GeV. When DM candidate can annihilate into two light Higgs bosons 
(mo > 120 GeV), one can derive the almost same A\ t D for three kinds of mixing cases and two 
R scenarios, which means that the light Higgs mixing and the parameter R do not significantly 
affect the total annihilation cross section. This conclusion can also be applied to 200 GeV < m D < 
500 GeV as shown in Figs. [7] and [F] (left panels). 

For the DM indirect search, the 2HBDM has two enhanced regions for (crv)o as shown in Fig. 
[6l Therefore the PAMELA antiproton measurements can exclude two very narrow regions. The 
predicted (ctv)q is the same as that in the 1HBDM for most parts of parameter space. When 



16 



200 GeV < m D < 500 GeV, one can still obtain <o-v> » 2.3 x 10~ 26 cm 3 sec" 1 . It is clear that 
different mixing cases and R scenarios lead to the same (crv)o except the resonance regions. 
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FIG. 6: The predicted thermally averaged DM annihilation cross section (o~v)o in the 2HBDM. 



In the 2HBDM, the DM-quark coupling a q in Eq. ([2TI) is given by 



M,pm q I f x f 3 f 5 
= -= Hr + -T + -T > (30) 



m 2 h m\ m\ 



2m D 

where /; have been presented in Appendix Eq. (IA3I) . Notice that we have neglected the parameters 
fi, Ia and f 6 since their contributions to erf 1 are velocity-dependent. Using the predicted A\^ D in 
Fig. [5](left panels), we calculate the spin-independent DM-nucleon elastic scattering cross section 
erf 1 for three mixing cases and two R scenarios. Different from (crv)o, the predicted a s n 1 obviously 
depends on the mixing and R as shown in Fig. [5] (right panels). Although three kinds of mixing 
cases have the almost same coupling A\ D for m D > 120 GeV in the R = 1 scenario, the predicted 
cr s n 1 in the Case III is far less than that in the Case I and Case II. This is because that there is 
cancellation between f\/m\ and f$lm\ in Eq. (1301) for the Case III. When the DM candidate can 
annihilate into two light Higgs bosons, a large R does not obviously affect the predicted coupling 
A^ D . However, the parameters /i,/3 and / 5 in Eq. (1301) will be significantly enlarged. Therefore 
erf 1 usually increases as R increases. The Case I clearly demonstrates this feature. The enlarged 
erf 1 in the R = 5 scenario may approach the CDMS II upper bound, which can be used to explain 



the two possible events observed by the CDMS II 112911 . It is worthwhile to stress that the Case II 
in the R = 5 scenario give a smaller erf 7 than that in the R = 1 scenario due to the cancellation 
from the different Higgs boson contributions. We conclude that the predicted erf 1 is significantly 
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dependent on the light Higgs mixing and the parameter R. For 200 GeV < m D < 500 GeV, the 
same conclusion can also be derived as shown in Figs. [7] and [8] (right panels). 

As shown in Figs. □ and [8] (right panels), the CDMS II U and XENON10 \m experiments 
can exclude the region m D < 50 GeV. For 200 GeV < m D < 500 GeV, our results show an upper 
bound for erf 1 which is still below the current experiment upper bounds. The future experiments 
XENON100 Q3JJ], CDMS 100 kg Q32D and XENON IT [33D can cover most parts of the allowed 
parameter space except the extreme cancellation cases. Nevertheless, it is still difficult to detect 
the DM direct or indirect signals for the resonance regions 50 GeV < m D < 60 GeV and 80 GeV < 
m D < 90 GeV. 



V. CONCLUSIONS 



In conclusion, we have investigated a scalar boson S u as the DM candidate in the left-right 
symmetric gauge model with two Higgs bidoublets, in which the SCPV can be easily realized. The 
stability of DM candidate S D is ensured by the fundamental symmetries P and CP of quantum field 
theory. In order to well understand the DM properties in the 2HBDM, we have firstly analyzed the 
1HBDM and shown that the predicted DM direct and indirect detection cross sections (erf/ and 
(crv)o) are the same as that in the minimal extension of SM with a real singlet scalar if m D < m h o. 
When the annihilation channel S D S D — » h°h° is open (m D > m h o), the H® exchange diagram 
relevant to A R leads to a continuous DM-nucleon elastic scattering cross sections erf/. Comparing 
with the 1HBDM, there are more scalar particles entering the DM annihilation and scattering 
processes in the 2HBDM. In the explicit calculations, we have considered three typical mixing 
cases and two Yukawa coupling scenarios (R = 1 and R = 5) to analyze the 2HBDM. It has 
been shown that (crv) is not sensitive to the light Higgs mixing and Yukawa couplings except the 
resonance regions. However erf 1 is significantly dependent on the above two factors. In general, 
erf 1 can be enhanced by large Yukawa couplings and approach the CDMS II upper bound, which 
can be used to explain the two possible events observed by CDMS II. It should be mentioned that a 
large Yukawa coupling may lead to a very small erf 1 in the extreme mixing case. Our results show 
that the future DM direct search experiments can cover most parts of the allowed parameter space. 
The PAMELA antiproton data can exclude two very narrow regions in the 2HBDM. In addition, 
we have shown that it is very difficult to detect the DM direct or indirect signals for the resonance 
regions since the Breit-Wigner resonance effect simultaneously suppresses erf 1 and (crv) . 
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Appendix A: Annihilation cross section 



For the annihilation processes So So — * ff, the annihilation cross section & f f is given by 



where 



with 



X 1 \ Am 2 



/ 



Pm = 



hi 



+ 



3,4 



+ 



ft 



5.6 



s - m 2 h + im h T h s - m 2 H + im H Y H s - m 2 A + im A Y A 



(Al) 



(A2) 



f\ = c x c z - Rc y s x - Rc x s y s z , f 2 = Rs x s y - Rc x c y s z , 
f 3 = Rc x c y + c z s x - Rs x s y s z , f 4 = -Rs x s z c y - Rc x s y , 
f 5 = Rs y c z + s z , f 6 = Rc y c z . 



(A3) 



The parameter R has been defined in Eq. (|29l) . The decay widths of three light neutral Higgs are 
given by 



h,H,A - 



&7TV 



h,H,A 



where jh = c 2 c 2 , y H = s 2 c 2 and y A = s 2 . ^ HA and T h l HA have the following forms: 

,3 



in 



h,H.A 



h,H,A 



327TV lw 
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Ami 
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h,H,A 



Ami 12«4 ^ 
l-^^ + ^i 
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h,H,A ^2 



EW 



4m?, 



h,H,A m h,H,A ) 



VV, 



III 



h,H,A 



Am 2 Wi ^ \2m A w ^ 



m 



h,H,A 



III 



h,H,A } 



For the annihilation processes S D S D — » Z\Z\ and S D S D — » W\ W\ , we have 



'l- 



Aml i 



X 



16^ \ 

2c x c z 



Ami \2m% \ £ 
s s 2 I A 
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2S X C 7 
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2s 7 



s - mf t + im h T h s - m 2 H + im H T H s -m 2 A + im A T A 
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4< ' 
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8tt > 

2c x c z 



+ 



Ami, 
1 ^ + 
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2s x c 7 



4 



+ 



2s, 



s - m 2 h + im h T h s - m 2 H + im H T H s -m 2 A + im A T A 



,(A4) 



(A5) 



(A6) 



(A7) 
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If the annihilation productions are a Higgs and a gauge boson, we can derive 



where 



cr w±h* 



1.5 



0"Z,A = 



1.5 



4z> [( j - m i- m l,) 2 - 4m >li] 

327T 5 
4fl [(* ~ m H ~ m l,) 2 ~ 4m ff m I,] 

32n s 

a Id [( s - m l- m V 2 - 4m l m2 z] 1 
32n s 

A\ D [{s - m 2 H± - m 2 Wi ) 2 - 4m 2 H± m 2 Wl ] 
An s 



2c 7 s, 
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2c z s x 
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(A8) 



Ct \ — CyC 7 XC £ Sy . 

@2 — C^ylCy $y) CyS x S-7 H~ is x SyS% , 

a 3 = ic y s x + s y s x - c x (c y s z - is y s z ) . 
When two DM candidates annihilate into two Higgs particles, we can obtain 
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o~kk 



cr,-; = 
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with 
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3ml 



mi 
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3mf 



+ 



m 



H 



3m\ 



m. 



<X\&XDV R 1 

+ ^- + 



S ~ m rjO h,D * - m l 

tl- 



(A9) 



(A10) 



(All) 



The subscripts k and ij run over (h, H, A) and (hH, hA, HA), respectively. The quantity F is defined 
as F(g) = arctanh(^)/^ with = ^jl - 4m 2 D /s ^J(s - mf - m 2 ) 2 - 4m 2 m 2 /(s - mf - m 2 ). The 
parameter fyj is given by y8,- 7 - = J(s - mf - m 2 ) 2 - 4mfm 2 / s. 
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